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Abstract 

The invariant-mass spectrum of e+e~ pairs produced in ^^C+^^C collisions at an incident en- 
ergy of 2 GeV per nucleon has been measured for the first time. The measured pair production 
probabilities span over five orders of magnitude from the vr'^-Dalitz to the pjio invariant-mass re- 
gion. Dalitz decays of -k^ and t] account for all the yield up to 0.15 GeV/c^, but for only about 
50% above this mass. A comparison with model calculations shows that the excess pair yield is 
likely due to baryon-resonance and vector-meson decays. Transport calculations based on vacuum 
spectral functions fail, however, to describe the entire mass region. 

PACS numbers: 25.75.-q, 25.75.Dw, 13.40.Hq 
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The properties of hot and dense hadronic matter represent a key problem in heavy- 
ion physics, with far-reaching imphcations for other fields such as the physics of compact 
stars. They are governed by non-perturbative QCD and cannot be derived directly from 
the underlying Lagrangian. Models predict, however, that hadron properties, such as mass 
and lifetime, depend on the temperature and density of the medium. While some hadronic 
many-body calculations give a broadening of the meson in-medium spectral function, other 
approaches predict dropping meson masses as precursors of chiral symmetry restoration l| . 

Experimentally, in-medium properties are difficult to observe. Suitable probes are dilep- 
tons or e^e~) from decays of short-lived resonances produced inside the hadronic 

matter created in the course of relativistic or ultra-relativistic heavy-ion collisions. At the 
CERN SPS, the CERES collaboration has established a significant excess of the dielectron 
yield as compared to that expected from the decays of hadrons after chemical freeze-out 
in 40 and 158 AGeV Pb+Au collisions j^. At 1 AGeV, the DLS collaboration found un- 
expectedly large electron-pair yields in C-l-C and Ca-|-Ca collisions [3]. In contrast to the 
situation at SPS beam energies, the DLS results cannot be described satisfactorily within the 
various scenarios proposed for possible changes of the in-medium spectral functions 
Indeed, the pair yields in the invariant-mass range between 0.15 and 0.6 GeV/c^, i.e. just 
below the p meson pole mass, still remain to be explained B, fl y, 3- Recently, the NA60 
collaboration has published new data on dimuon production which allow a determination of 
the in-medium spectral function of the p meson in 158 AGeV In+In collisions lOj. 

The High-Acceptance DiElectron Spectrometer HADES at GSI, Darmstadt, operates in 
the SIS/Bevalac energy regime of 1-2 AGeV. In this Letter we report on the first measure- 
ment of inclusive electron-pair production in ^^C+^^C collisions at a kinetic beam energy of 
2 AGeV. 

A carbon beam of 10^ particles/s was incident on a 2-fold segmented carbon target 
with a thickness corresponding to 2 x 2.5% interaction lengths. The HADES spectrometer. 



described in detail in Refs. 
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121], consists of a 6-coil toroidal magnet centered on the 



beam axis and six identical detection sections located between the coils and covering polar 
angles between 18° and 85°. In the measurement presented here, each sector was composed 
of a gaseous Ring-Imaging Cherenkov (RICH) detector, two planes of Mini-Drift Chambers 
(MDC-I and MDC-II) for track reconstruction and a Time-Of-Flight wall (TOF/TOFino) 
supplemented at forward polar angles with Pre-SHOWER chambers. The interaction time 
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was obtained from a fast diamond start detector located upstream of the target. This 
geometry results in a smooth dielectron acceptance, shown in Fig. [1] as a function of 
invariant mass and transverse momentum, averaged over the rapidity range < y < 2. 
The data readout was started by a first-level trigger (LVLl) decision, requiring a charged- 
particle multiplicity MUL > 4 in the TOF/TOFINO detectors, accepting 60% of the total 
cross section. It was followed by a second-level trigger (LVL2) requesting at least one electron 
track. With this trigger condition, a 10-fold pair enrichment at a pair efficiency of > 92% 
was achieved. Furthermore, the LVL2 introduced no bias on the shapes of measured pair 
distributions, as checked by a direct comparison to the unbiased LVLl events. 

The results presented below were obtained from events with a positive LVL2 decision with 
a total statistics corresponding to 6.5 x 10^ LVLl events. The electron track reconstruction 
proceeded in four steps (see also |l2l]): 

(1) The RICH ensures electron detection for momenta p > 0.05 GeV/c. From the mea- 
sured ring positions, polar (9) and azimuthal (0) angles are calculated assuming electron 
emission from the target. 

(2) The 6, values of electron candidates in the RICH were correlated with the track 
angles reconstructed in the MDC within windows corresponding to ±2 standard deviations, 
as derived from the electron distributions. 

(3) Candidate tracks were then matched with those hits in the TOF or TOFINO/Pre- 
SHOWER detectors fulfilling electron conditions, i.e. (a) a particle velocity of /3 = 1 ± 3cr/3, 
with a/3 given by the time-of-fiight resolution, and (b) an electromagnetic shower signal in 
the Pre-SHOWER. 

(4) Finally, the track momentum was determined by a fit of an appropriate track model 
to the reconstructed hit positions making use of the defiection in the known magnetic field. 

The identified single-electron tracks were combined into opposite-sign pairs from which 
invariant-mass distributions, with the mass resolution (JMee/^ee = 9% at M^e = 0.8 GeV/c^, 
were built. Many of these e+e~ pairs, however, represent combinatorial background (CB) 
which has to be reduced and subtracted. The CB is mostly due to uncorrelated electrons 
from TT^ —>■ 77 decays followed by photon conversion, either in the target or in the RICH 
radiator, and/or from vr'' e^e~7 Dalitz decays. Such pairs have small opening angles 
and often produce partially overlapping tracks in the MDC. They are rejected efficiently 
by applying conditions on the opening angle, 9e+e- > 9°, and on the fit quality (x^) of 



5 




FIG. 1: Geometrical acceptance of the HADES spectrometer for e~^e pairs with laboratory opening 
angle O^+e- > 9*^ as a function of their invariant mass and transverse momentum. 

the reconstructed track segments, removing 95% of the conversion pairs while reducing the 
dielectron signal with Mee > 0.15 GeV/c^ by less than 10%. 

Fig. [2] shows the resulting e'^e" invariant-mass distribution, its signal part, the CB and 
the resulting signal/CB ratio. Like-sign e'^e'^ and e~e" pairs were formed and subjected 
to the same selection criteria as the opposite-sign pairs. From the reconstructed like-sign 
invariant-mass distributions, dN^^ / dM^e and dN /dMee-, the respective CB distribution 
was calculated as Nqb = 2-\/ N^^N . For masses Mee > 0.5 GeV/c^, where statistics is 
smaller, the CB was obtained by an event mixing procedure. Uncorrelated opposite-sign 
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e'^e~ pairs were formed from different events but originating from reactions in the same 
target segment. It was verified that for Mge > 0.15 GeV/c^ the CB distributions obtained 
from the event mixing and the hke-sign pairs agree within 10%. All distributions were 
normalized to the number of neutral pions A^^o (see below). In total ~ 23000 signal pairs 
(~ 2000 with Mee > 0.15 GeV/c^) were reconstructed. 
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FIG. 2: Reconstructed (not efficiency corrected) e~^e~ invariant-mass distributions (open circles) for 
pairs with opening angles ^?e+e- > 9°. The signal distribution (dots) was obtained by subtracting 
the combinatorial background (full triangles). Errors indicated are statistical only. Inset: Signal- 
to-background (S/B) ratio versus invariant mass. 



We corrected the spectra for detector and reconstruction inefficiencies by Monte-Carlo 



simulations embedding electron tracks with uniform 1/p and isotropic angular distributions 
into ^^C+^^C events generated with the UrQMD transport model [l^. Resulting events 
were digitized and processed through the same analysis chain as the measured data. The 
single-electron efficiencies, e±, were calculated as a function of charge (±), momentum (p), 
polar {9) and azimuthal (0) emission angles. The data were corrected on a pair-by-pair 
basis with the weighting factor with = e+ ■ e_ for given electron momenta and 

emission angles. The CB was treated likewise and subtracted, as described above, to obtain 
the efficiency-corrected pair signal distribution. This prescription relies on the assumption 
that the single- leg efficiencies are independent, as was carefully checked in our simulations 
and proven to be valid within 15% for pairs with opening angles ^e+e- > 9*^. The geometrical 
pair acceptance of the HADES detector was obtained in analogy to the pair efficiency as 
the product of two single-electron acceptances A±{j),9,(f)). The resulting matrices, together 
with a momentum resolution function, constitute the HADES acceptance filter (available 
upon request). We made no attempt to extrapolate the measured dielectron yields to the 
full solid angle. 

Fig. [3^ shows the e^e~ invariant-mass distribution of the signal pairs after efficiency 
correction and normalized to the average number of charged pions N^^ = ^{Nt^+ + N^^-). 
The latter were identified in HADES by means of the time-of-fiight measurement [l^ and 
their yield was extrapolated to full solid angle, taking into account our measured angular 
distributions found to be in agreement with UrQMD calculations. In the isospin-symmetric 
system ^^C-|-^^C, is in fact a good measure of the tt" yield. This way of normalizing the 
pair spectra compensates to ffist order the bias caused by the implicit centrality selection 
of our trigger. Indeed, simulations based on UrQMD events indicate that LVL2 events have 
an average number of participating nucleons Apart = 9.0, instead of 6 for true minimum-bias 
events. The pion multiplicity per number of participating nucleons M,,/ Apart = 0.137±0.015 
obtained in our experiment agrees with previous measurements of charged and neutral pions 



151] within the quoted error of 11%. The error is dominated by systematic uncertainties 



in the acceptance and efficiency corrections of the charged-pion analysis and represents our 
overall normalization error. 

We ffist compare our results with a pair cocktail (cocktail A) calculated from free 7r°, 
7] and uj meson decays only (Fig[3K). This cocktail aims at representing all contributions 
emitted after the chemical freeze-out of the ffieball. While the ffist two sources are directly 
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ansatz 



constrained by data [IJ] , the production rate of the uj meson is taken from an m^-scahng 



161 ] ■ In our event generator (PLUTO 17|]) meson production was hence modeled 



assuming emission from a thermal source with a temperature T = 80 MeV, but no radial 
expansion velocity (/3,.=0). Furthermore, for the vr" mesons, an anisotropic angular distri- 
bution of the type dN/dcos{9cM) ~ 1 + 02 ■ cos'^{9cm) with 02 = 0.7 was used, as deduced 
from our charged pion analysis. The accepted lepton-pair yield was found to change by less 
than 14% when varying the source parameters over a broad range: jSr (0-0.3), 02 (0-1.0) and 
T (50-90 MeV). 

While experimental data and simulated cocktail A (solid line in Fig. [3K) are in good 
agreement in the tt" region, the coctail undershoots the data for M > 0.15 GeV/c^ and 
clearly calls for additional sources. Such contributions are indeed expected from the decay 
of short-lived resonances, mainly the A(1232) and the p, excited in the early phase of the 
collision. To include in our cocktail pairs from A'''^ Ne^e" decays, we assumed that 



the A yield scales with the vr yield and employed a calculated decay rate .9!]. To add the 
p meson we used a similar prescription as for the uj. For this broad resonance {T p = 150 
MeV), m_L scaling strongly enhances the low- mass tail, resulting in the skewed shape visible 
in Fig. [3^. The long-dashed line in this figure shows the comparison of the full thermal 
cocktail (cocktail B) with our data. As expected, the simulated yield above 0.15 GeV/c^ is 
now increased and, in particular, the high-mass region is populated with dielectrons from 
p ^ e^e" decays, but the calculation still falls short of reproducing the data. 

To discuss in more detail the excess pair yield, we show in Fig. [3b the ratio of the data and 
cocktail A. Statistical and systematical errors originating from the CB subtraction, efficiency 
corrections and normalization, are shown separately. In the intermediate mass range of 
0.15 — 0.50 GeV/c^, the enhancement factor above the dominant rj contribution amounts to 
F(2.0) = Ftot(2.0)/r^(2.0) = 2.07 ± 0.21{stat) ± 0.38{sys). Assuming that the excess pairs 
have in this mass region an overall acceptance close to the one of eta Dalitz pairs, one can 
compare F[2.0) to the enhancement factor measured in C+C by DLS at a beam energy of 
1.04 AGeV 3[. Using the DLS data and an appropriately filtered PLUTO cocktail generated 
for this bombarding energy, we obtain a factor of F(1.04) = 6.5 ± 0.5{stat) ± 2.1{sys). 
Knowing that, going from 1.04 to 2 AGeV, inclusive 77 production in C-l-C collisions increases 
by a factor y^(2.0)/F^(1.04) = 13 ± 3 [ij, [l8|, the energy scaling factor of the excess pair 
yield y'ea;c(2.0)/Yesc(l-04), where Fexc = ^tot — ^r?, can be deduced from the two enhancement 
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factors F(2.0) and F(1.04). It follows that Ye^c{'2.0)/Ye^c{lM) = 2.5±0.5(stot) ± 1.5(s?/s). 
This energy scaling is remarkably similar to the known scaling of pion production, i.e. 
Y„{2.0)/Y^{1.0A) = 2.3 ± 0.3 Q, [l^. It suggests that the pair excess is indeed driven by 
pion dynamics, involving e.g. A and p excitations. 

At still higher masses the ratio of data and cocktail A develops a pronounced maximum 
around Mge ~ 0.6 GeV/c^, mainly due to the lack of p decays in cocktail A. Especially the 
mass region between the 1] and the u pole is expected to be dominated by the thermally 
populated low- mass tail of the broad p resonance. Cocktail B (long-dashed line in Fig. [3)d), 
which includes p and A decays, shows indeed an enhancement, but still does not fully explain 
the observed pair yield. 

Beyond the cocktail calculations discussed above we have also made a comparison with 
various transport models. The latter, besides offering a realistic treatment of collisions 
dynamics, also handle the propagation of broad resonances, related off-shell effects and 
multi-step processes, all known to play a crucial role at our bombard.ing energy. Hence, 
dielectron distributions were calculated with the HSD [6], RQMD [8] and UrQMD [3] 
transport models (assuming vacuum spectral functions). All calculations were filtered with 
the HADES acceptance and normalized to their respective 7r° yields. The respective ratios 
of calculated dielectron yields to our cocktail A are shown as curves in Fig. [3)d. Due to 
the normalization, all models agree in the 7r° mass region and give ratios consistent with 
unity. In the intermediate mass region all transport models fall short, with UrQMD coming 
closest to the data. As checked, all models agree on the r] contribution itself within 20%, 
and the discrepancy can be traced to differences in their respective treatment of population 
and decay of the baryonic resonances [g, y, Q. At higher masses, M^e > 0.5 GeV/c^, all 
models qualitatively reproduce the trend given by the data, but differ in the yield due to 
different amplitudes used for the couplings to intermediate resonances j^, Q • All transport 
models consistently overestimate the pair yield at the p, u meson poles. Final conclusions 
can be drawn only when more refined calculations with in-medium spectral functions become 
available. 

In summary, we report on the first measurement of inclusive dielectron production in 
i2q_|_i2q collisions at £'^^^^=2 AGeV, spanning five orders of magnitude in yield. At low 
masses, i.e. Mee < 0.15 GeV/c^, the pair yield is in agreement with the known vr" production 
and decay probabilities. For 0.15 GeV/c^ < Mge < 0.5 GeV/c^ it exceeds expectations based 
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on the known production and decay rates of the f] meson by 2.07±0.21(stat)±0.38(st/s). This 
pair yield excess is consistent with that measured by DLS at 1.04 GeV if its energy scahng 
is similar to that of pion production. Additional sources associated with the radiation from 
the early collision phase (A°(+^ Ne~^e~ , p e^e~) are needed to account for the excess 
observed for M > 0.15 GeV/c^. However, transport calculations based on vacuum spectral 
functions only still fail to quantitatively describe the excess yield in the full invariant-mass 
range. 
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FIG. 3: (a) Efficiency- and background-corrected e^e^ invariant-mass distribution for Oe+e- > 
9" (symbols) compared to a thermal dielectron cocktail of free vr'^, rj and cu decays (cocktail A, 
solid line), as well as including p and A resonance decays (cocktail B, long-dashed line). Only 
statistical errors are shown, (b) Ratio of data and cocktail A (dots), compared to ratios of various 
model calculations and cocktail A. All calculations have been filtered and folded with the HADES 
acceptance and mass resolution. Statistical and systematic errors of the measurement are shown 
as vertical and horizontal lines, respectively. The overall normalization error of 11% is depicted by 
the shaded area. 13 



